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O ! 

, We study an exactly solvable version of the famous random Boolean satisfiability problem, the 

^N) ' so called random XOR-SAT problem. Rare events are shown to affect the combinatorial "phase 

; I ' diagram" leading to a coexistence of solvable and unsolvable instances of the combinatorial problem 

, in a certain region of the parameters characterizing the model. Such instances differ by a non- 

extensive quantity in the ground state energy of the associated diluted spin-glass model. We also 
show that the critical exponent v, controlling the size of the critical window where the probability 
lO ■ of having solutions vanishes, depends on the model parameters, shedding light on the link between 

random hyper- graph topology and universality classes. In the case of random satisfiability, a similar 
behavior was conjectured to be connected to the onset of computational intractability. 



C/3 



o 
o 



X 



PACS Numbers : 89.80.-fh, 75.10.Nr 



I. INTRODUCTION 



The satisfaction of constrained Boolean formulae, the so called Satisfiability (SAT) problem, is a key problem of 
_JL ' complexity theory in computer science that can be recast as an energy minimization problem (ground state search) 
^ ' in diluted spin glass models. Many hard computational problems have been shown to be NP-complete Q through a 
Q ' polynomial mapping onto the SAT problem, which in turn was the first problem identified as NP-complete by Cook 
O ; in 1971 §. 

. Recently Q , the research activity has become more and more focused on the study of the random version of SAT 
problem defined as follows. Consider TV Boolean variables Xi, i = l,...,iV. Call clause C the logical OR of K 
^ , randomly chosen variables, each of them being negated or left unchanged with equal probabilities. Then repeat this 
' process by drawing independently M random clauses C^, a — 1, . . . , M. The logical AND of all clauses is said to 
' be satisfiable if there exists a logical assignment of the {xi} evaluating ^ to true, unsatisfiable otherwise. 
■ Numerical experiments have concentrated upon the study of the probability Pn{j,K) that a given J- including 
[ M = clauses be satisfiable. For large sizes, there appears a remarkable behavior : Poo{"/,K) seems to be unity 
. for 7 < "fc{K) and vanishes for 7 > "fdK) Such an abrupt threshold behavior, separating a so-called SAT phase 

from an UNSAT one, has indeed been rigorously confirmed for 2-SAT, which is in P, with 7c(2) = 1 For K > 3, 
"p-, K-SAT is NP-complete and much less is known. The existence of a sharp transition has not been rigorously proven 
yet but relatively good estimates of the thresholds have been found: 7c (3) ~ 4.2 4.3 0. Moreover, some rigorous 
^ , lower and upper bounds to 7c (3), have been established p. 

I ' The interest in random K-SAT arises from the fact that it has been observed numerically that hard random instances 
I are created when the problems are critically constrained, i.e. close to the SAT/UNSAT phase boundary [||J|]. The 
^ i study of such hard instances represent a theoretical challenge towards a concrete understanding of complexity and 
the analysis of algorithms. Moreover, hard random instances are also test-bed for the optimization of heuristic 
(incomplete) search procedures which are widely used in practice. 

The statistical mechanics study of random K-SAT have provided some geometrical understanding of the onset of 
complexity at the phase transition through the introduction of a functional order parameter which describes the 
5^ ■ geometrical structure of the space of solutions. The nature of the SAT/UNSAT transition for the different values 
\ of K appears to be a particularly relevant prediction Q. The SAT/UNSAT transition is accompanied by a smooth 
(respectively abrupt) change in the structure of the solutions of the 2-SAT (resp. 3-SAT) problem. More specifically, 
at the phase boundary a finite fraction of the variables become fully constrained while the entropy density remains 
finite. Such a fraction of frozen variables (i.e. those variables which take the same value in all solutions) may undergo 
a continuous (2-SAT) or discontinuous (3-SAT) growth at the critical point. This discrepancy is responsible for the 
difference of typical complexities of both models recently observed in numerical studies. The typical solving time of 
search algorithms displays an easy-hard pattern as a function of 7 with a peak of complexity close to the threshold. 
The peak in search cost seems to scale polynomially with N for the 2-SAT problem and exponentially with N in the 
3-SAT case. From an intuitive point of view, the search for solutions ought to be more time-consuming in presence of 
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a finite fraction of fully quenched variables since the exact determination of the latter requires an almost exhaustive 
enumeration of their configurations. 

To test this conjecture, a mixed 2 + p-model has been proposed, including a fraction p (resp. 1 — p) of clauses 
of length two (resp. three) and thus interpolating between the 2-SAT {p — 0) and 3-SAT {p — 1) problems. The 
statistical mechanics analysis predicts that the SAT/UNSAT transition becomes abrupt when p > po — 0.4 [p|-pT|. 
Precise numerical simulations support the conjecture that the polynomial/exponential crossover occurs at the same 
critical po- Though the problem is both critical (7c = 1/(1 — p) for p < po) and NP-complete for any p > 0, it is 
only when the phase transition becomes of the same type of the 3-SAT case that hardness shows up. An additional 
argument in favor of this conclusion is given by the analysis of the finite-size effects on ^^(7, K) and the emergence of 
some universality for p < pq. A detailed account of these findings may be found in |-|l|. For p < Pq the exponent v, 
which describes the shrinking of the critical window where the transition takes place, is observed to remain constant 
and close to the value expected for 2-SAT. The critical behavior is the same of the percolation transition in random 
graphs (see also ref. [p^). For p > po the size of the window shrinks following some p-dependent exponents toward 
its statistical lower bound but numerical data did not allow for any precise estimate. 

In this paper, we study an exactly solvable version of the random 2-|-p SAT model which displays new features and 
allows us to settle the issue of universality of the critical exponents. The threshold of the model can be computed 
exactly as a function of the mixing parameter p in the whole range p G [0, 1]. Rare events are found to be dominant 
also in the low 7 phase, where a coexistence of satisfiable and unsatisfiable instances is found. A detailed analysis for 
the p — 1 case can be found in ref. 15|. 

The existence of a global - polynomial time - algorithm for determining satisfiability allows us to perform a finite 
size scaling analysis around the exactly known critical points over huge samples and to show that indeed the exponent 
controlling the size of the critical window ceases to maintain its constant value 1/ = 3 and becomes dependent on p as 
soon as the phase transition becomes discontinuous, i.e. for p > po = .25. Above po and below pi ~ 0.5, the exponent 
v takes intermediate values between 3 and 2. Finally, above pi the critical window is determined by the statistical 
fluctuations of the quenched disorder H and so 2. 



II. MODEL DEFINITION AND OUTLINE OF SOME RESULTS 

The model we study can be viewed as the mixed 2 +p extension of the 3-hyper-SAT (hSAT) model discussed in , 
as much as the 2-(-p-SAT |^ is an extension of the usual K-SAT model. In computer science literature the hSAT model 
is also named XOR-SAT and its critical behavior is considered an open issue ||l^. Given a set of N Boolean variables 
{xi = 0, l}i=i,....Ar we can write an instance of our model as follows. Firstly we define the elementary constraints (a 
mixture of 4 and 2-clauses sets with 50% satisfying assignments): 

C{ijk\ + 1) = {xi V Xj V Xk) A {xi V Xj V Xk) A {xi V Xj V Xk) A {xi V Xj V Xk) 

C{ijk\ — 1) = {xi V Xj V Xk) A {xi V Xj V Xk) A {xi V Xj V Xk) A {xi V Xj V Xk) , (1) 
for the 3-hSAT part, and 

C{ij\ + 1) ^ {xi V Xj) A (ij V Xj) 

C{ij\ - 1) = {x, V Xj) A (5, V Xj) , (2) 

for the 2-hSAT part. A and V are the logical AND and OR operations respectively and the over-bar is the logical 
negation. A more compact definition can be achieved by the use of the exclusive OR operator ©, e.g. C{ijk\ + I) = 
Xi (B Xj (B Xk- Then, we randomly choose two independent sets E3 and E2 of pM triples k} and (1 — p)M couples 
{i,j} among the N possible indices and respectively pM and (1 —p)M associated unbiased and independent random 
variables Tijk = ±1 and Jij = ±1, and we construct a Boolean expression in Conjunctive Normal Form (CNF) as 

F= /\ C{ijk\T,jk) A ^feV^j) • (3) 
{ij,fc}eB3 {«j}e-E2 

As in ||l5|, we can build a satisfiable version of the model choosing clauses only of the C{ij\ + 1) and C{ijk\ + 1) type. 
For p < Po the problem is easily solved by local and global algorithms, whereas interesting behaviors are found for 
p > Po, where the local algorithms fail. 

The above combinatorial definition can be recast in a simpler form as a minimization problem of a cost-energy 
function on a topological structure which is a mixture of a random graph (2-spin links) and hyper-graph (3-spin 
hyper-links). We end up with a diluted spin model where the Hamiltonian reads 
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{i,].k}£E3 {i.]}eE2 

where the 5*^ are binary spin variables and the the random couphngs can be either ±1 at random. The satisfiable 
version is nothing but the ferromagnetic model: T^tj = 1 and Jy — 1 for any link. 

As the average connectivity 7 of the underlying mixed graph grows beyond a critical value 7c (p), the frustrated 
model undergoes a phase transition from a mixed phase in which satisfiable instances and unsatisfiable ones coexist to 
a phase in which all instances are unsatisfiable. At the same 7c(p) the associated spin glass system, undergoes a zero 
temperature glass transition where frustration becomes effective and the ground state energy is no longer the lowest 
one (i.e. that with all the interactions satisfied). At the same critical point the unfrustrated, i.e. ferromagnetic, version 
undergoes a para-ferro transition, because the same topological constraints that drive the glass (mixed SAT/UNSAT 
to UNSAT) transition in the frustrated model are shown to be the ones responsible for the appearance of a nonzero 
value of the magnetization in the unfrustrated one . We shall take advantage of such coincidence of critical lines 
by making the analytical calculation for the simpler ferromagnetic model. 

Moreover, the nature of the phase transition changes from second to random first order, when p crosses the critical 
value Pq = 1/4. For p > pq the critical point 7^ {p) is preceded by a dynamical glass transition at 7^ {p) where ergodicity 
breaks down and local algorithms get stuck (local algorithms are procedures which update the system configuration 
only by changing a finite number of variable at the same time, e.g. all single or multi spin flip dynamics, together 
with usual computer scientists heuristic algorithms). The dynamical glass transition exist for both versions of the 
model and corresponds to the formation of a locally stable ferromagnetic solution in the unfrustrated model 
(the local stability is intimately related to the ergodicity breaking). 

III. STATISTICAL MECHANICS ANALYSIS 

Following the approach of ref. we compute the free energy of the model with the replica method, exploiting the 
identity log ^ Z" ':^= 1 + n <C logZ ^ +0(n^). The n*^ moment of the partition function is obtained by replicating 
n times the sum over the spin configurations and then averaging over the quenched disorder 

Y «exp [ -/3^ffj[S''] ] » , (5) 

Si,S^...,S" \ a=l / 

Since each of the M clauses is independent, the probability distributions of the ferromagnetic couplings can be written 
as 



pmjk})^ n 



i<j<k 



giving the following expression for the <C Z'"- ^ 



P<„,„.n[(i-^)«,)+^ 



(6) 



«Z"»= Y exp|-^7iVn-7A^+f^E^'^"'''^"'^" + ^^E^'^"'''^"+0(l)^ (7) 
Introducing the occupation fractions c(<t) (fraction of sites with replica vector a) , one gets 

- /3F[c] = -7(1 + M - E c{a) log c(a) + {l-phJ2 """^ + PI E c(a)c(p)c(f)e^ 2. ^"""^^ . (8) 

In the thermodynamic hmit we can calculate the free energy via the saddle point equation obtaining 

c{a) = exp J -A + 2(1 - p)iY.<P) ^MPY.'^'p") + 3p7E exp(/3 E '^V"^") \ (9) 
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The Lagrange multiplier A = —7(2 + p) ensures the normalization constraint J^ff'^i'^) = 1 in the limit n ^ 0. 
Finding the minimal (zero in the unfrustrated case) value of the cost function amounts to studying the /3 — > 00 (zero 
temperature) properties of the model. In the Replica Symmetric (RS) Ansatz, the behavior of the spin magnetization 
can be described in terms of effective fields m = tanh f3h whose probability distribution is defined through 

-(")-y.j^^(^) (2cosh(/?/.))» ■ (^0) 

In the unfrustrated or ferromagnetic case, the P{h) turns out to have the following simple form 

P{h)=J2riS{h-l) , (11) 

where the effective fields only assume integer values. In the satisfiahle model the saddle point equations all collapse 
in one single self-consistency equation for tq: 

^ g-3p7(i-'-o)"-2(i-p)7(i-'-o) = ^ ^ ^-3p7^-2(i-p)7 (3p7)i (2(1 - ^7)2 _ (1 _ ^^^f^ [rof^ . (12) 

Ci=OC2=0 

The equations for the frequency weights ri with I > follow from the one for tq and read 

_ [3p7(l->-o)^ + 2(l-p)7(l-ro)]' ^ ^^^^ 

The previous self consistency equations for ro (or for the magnetization to = 1 — ro) can easily be derived by the 
same probabilistic argument used in Jisf , due to the fact that the clause independence allows to treat the graph and 
the hyper-graph part separately. Note that in the simple limit p — Q we retrieve the equation for the percolation 
threshold in a random graph of connectivity 7 JTot , 

00 \^ 

l-r, = e-'^Y.^^i^--o) ■ (14) 

k=0 

Since the ground state energy of the ferromagnetic model is zero, the free energy coincides with the ground state 
entropy, which can be written as a function of p, r^ and 7: 

^(7) = log(2)[ro(l - log(ro)) - 7(1 - p)(l - (1 - rof) ~ 7P(1 - (1 - ^0)')] (15) 

To find the value of the paramagnetic entropy we put ourself in the phase where all sets of 4- and 2-clauses act 
independently, each therefore dividing the number of allowed variables choice by two: the number of ground states 
wiU be Ng, = 2f^-P'rN-{i-phN ^ 2^(i-'>'). The resulting value of Spara = (1 - 7) log(2) coincides with the one found 
setting To = 1 in Eq.(^5|). This may not be the case in more complicated models, where the ground state entropy is 
a complicated function of 7 also for 7 < 7c, reflecting the fact that the magnetization probability distribution in the 
paramagnetic phase could be different from a single delta peak in m = 0. 
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Solving the saddle point equation for tq, we find that a paramagnetic solution with ro = 1 always exists, while at 
a value of 7 = 7d(p) there appears a ferromagnetic solution in the satisfiable model. For p = 0, the critical value 
coincides as expected with the percolation threshold jd{0) = 1/2. As long as the model remains like 2-SAT, up to 
p < Po = 0.25, the threshold is the point where the ferromagnetic solution appears and also where its entropy exceeds 
the paramagnetic one. The critical magnetization is zero and the transition is continuous. For larger values of the 
control parameter p the transition becomes discontinuous. There appears a dynamical transition at 7 = ^d{p) where 
locally stable solutions appear. At 7 = 7c (p) > ld{p), the non trivial ro 7^ 1 solution acquires an entropy larger than 
the paramagnetic one and becomes globally stable. The shape of 7 = 7d(p) and 7 — ^c{p) as functions of p are shown 
in Fig. 1^. The inset picture shows the magnetization of the model at the points where the dynamical and the static 
transitions take place. 



IV. NUMERICAL SIMULATIONS 



The model can be efficiently solved by a polynomial algorithm based on a representation modulo two (i.e. in Galois 
field GF[2]). If a formula can be satisfied, then a solution to the following set of M equations in variables exists 

S^SjSk = Tyfc y{i,j,k}eE^ , , 

S^SJ = J,, V{^,J}e£;2 ^""^ 

Through the mapping S^ = {-1)"', = (-1)'''^ and T^tc = with a^,r|^■jk,C^]k G {0,1}, Eq.@ can be 

rewritten as a set of binary linear equations 



(fTi + (7j + ak) mod 2 = djk V{i, j, k} £ E3 
{ai + aj) mod 2 = r]ij V{«, j} £ E2 



(17) 



For any given set of couplings {r]ij,Cijk}, the solutions to these equations can be easily found in polynomial time by 
e.g. Gaussian substitution. The solution to the M linear equations in A^ variables can be summarized as follows: a 
number Ndep of variables is completely determined by the values of the coupling {rjij, Cijfc} ^nd by the values of the 
Nfree = N — Ndep independent variables. The number of solutions is 2^J'''" and the entropy S{'^) — \og{2)N free/N = 
log(2)(l — Ndep{'y)/N). As long as Ndep = M we have the paramagnetic entropy Spara = log(2)(l — 7). However Ndep 
may be less than M when the interactions are such that one can generate linear combinations of equations where no 
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cr's appear, like = /{{rjijXijk})- This kind of equations correspond to the presence of loops (resp. hyper-loops [p^) 
in the underlying graph (resp. hyper-graph). A hyper- loops (generalization of a loop on a hyper-graph) is defined as 
a set S of (hyper-)links such that every spin (i.e. node) is "touched" by an even number of (hyper-)links belonging to 
S (see Fig. I). 




hyper-loop 
p>0 



FIG. 2. Typical loop and hyper-loop. Lines are 2-spin links, while triangles are 3-spin links. Note that every vertex has an 
even degree. 



Here we are interested in the fraction of satisfiable instances Psat{1tP)i averaged over the random couplings 
distribution. One can show that, for any random (hyper-)graph, Psat is given by 2~^'^', where Nhi is the number 
of independent (hyper-)loops In Fig. ^ we show the fraction of satisfiable instances as a function of 7 for 

p ~ and p = 0.5. The vertical lines report the analytical predictions for the critical points, jdp = 0) = 0.5 and 
7(p = 0.5) = 0.810343. 




I ' ' ^ ^-^-^ ' I - 1 

0.2 0.4 0.6 0.8 1 

Y 

FIG. 3. SAT probabilities PsAT{'y,p) for p = and p — 0.5. Data has been averaged over 10* different random hyper-graphs. 
Vertical straight lines are analytical predictions for critical points: jdp = 0) = 0.5 and 7c(p ~ 0.5) — 0.810343. Bold curves 
for 7 < 7c are analytical predictions for the SAT probability in the large A'^ limit. 



In the limit of large A'^ and for p — 0.5 the fraction of SAT instances sharply vanishes at the critical point in a 
discontinuous way, that is lim^_^^- PsAril) > while lim^^^+ Psat{i) = 0. This is the usual behavior already 
measured in 3-SAT |^,^ and 3-hyper-SAT [^, with the SAT probabilities measured on finite systems crossing at jc 
and becoming sharper and sharper as N increases. On the contrary for p — and large A^ the probability of being 
SAT becomes zero at 7c in a continuous way. The main consequence is that finite size corrections make -Psat (7) 
larger than its thermodynamical limit both before and after the critical point and thus the data crossing is completely 
missing. 

Note also that for p < 1 the fraction of SAT instances for 7 < jdp) is finite and less than 1 even in the thermo- 
dynamical limit, implying a mixed phase of SAT and UNSAT instances. This is due to the presence in the random 
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hyper-graph of loops made only by 2-spin links (indeed the mixed phase is absent for p = 1 when only 3-spin interac- 
tions are allowed p5| ). The expression for the SAT probability in the thermodynamical limit (bold curves in Fig. 
the lower most for p = and the uppermost for p = 0.5) can be calculated analytically and the final result is 

^SAT(7,P) = e^^(i-P)[i+^(i-f)l [1-27(1 for 7<7c(p) • (18) 

In order to obtain to above expression we note that the SAT probabihty is related to the number of (hyper-)loops by 

oc 

Psat(7,p)= E^("^;^'P)2"" ' (19) 

m=0 

where P(rn;j,p) is the probability of having m (hyper-)loops in a random (hyper-)graph with parameters 7 and p, 
and the factor 2~™ comes from the probability that for all the m (hyper-)loops the product of the interactions is 1 
(thus giving no contradiction in the formula). In order to estimate P(rn;7,p) we may restrict ourselves to consider 
only simple loops (made of 2-spin links), because hyper- loops which involve at least one 3-spin link are irrelevant in 
the thermodynamical limit. This can be easily understood with the help of the following counting argument. 

The probability that a given 2-spin link is present in a random (7,^) hyper-graph is p2 = 27(1 — p)/N and for a 
3-spin hyper-link is ps — 6jp/N'^. Thus the probability of finding in a random (7,p) hyper-graph a hyper-loop made 
of n2 links and ris hyper-links (na must be even) is just the number of different ways one can choose the (hyper-)links 
times P2^p^^ ■ Because the number of nodes belonging to a hyper-loop of this kind is at most n„ = 71,2 + 3^3/2 and 
the number of different hyper- loops of this kind is order iV"" , we have that the probability of having a hyper- loop 
with 712 links and hyper-links is order iV~"^/^. 

Then, for 7 < 7c the number of hyper-loops is still finite (their number becomes infinite only at 7c where a transition 
to a completely UNSAT phase takes place) and the SAT probability, in the large N limit, is completely determined 
by the number of simple loops (71.3 = 0). 

The typical number of these loops does not vanish for 7 < 2{i-p) ' therefore such "rare" events lead to a 
coexistence of SAT and UNSAT instances with equal energy density. 

The average number of loops of length k can be easily calculated and it is given by x'' / {2k), where x — 27(1 — p). 
The average number of loops of any size 

oc ^ 2 

k=3 

indeed diverges for a; ^ 1, that is for 7 2(i-p) ■ "^^^ probability of having m loops in a random (7,p) hyper-graph 
is then 



P(m;7,p) = e-^(-)^ 



(21) 



and the fraction of SAT instances turns out to be the one in Eq.(|l8|). 

We have numerically calculated the SAT probabilities for many p and N values, finding a transition from a mixed 
to a completely UNSAT phase at the 7c (p) analytically calculated in the previous section. We also find, in agreement 
with analytical results, that the transition is continuous as long as p < 1/4 and then it becomes discontinuous in the 
SAT probability. 

Let us now concentrate on the scaling with N of the critical region. We have considered several alternative definitions 
for the critical region. The one we present here seems to be the simplest and also the most robust, in the sense it can 
be safely used when the transition is both continuous (p < 0.25) and discontinuous {p > 0.25). We assume that the 
size of the critical region is inversely proportional to the derivative of the SAT probability at the critical point 



u^iN,p)-' = aP5AT(7,p) 



(22) 



7=7c 

For any value of p the width w{N) goes to zero for large N and the scaling exponent i>{p) is defined through 

wiN,p) o^N-^/^^P^ . (23) 
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FIG. 4. Scaling of the critical window width. Errors are smaller than symbols. Lines are fits to the data. 
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FIG. 5. Critical u exponents obtained from the fits shown in Fig. ^ For p — 0.75 and p - 
term power exponent, the leading term one being fixed to —1/2 (filled circles). 



I filled squares show the subleading 



In Fig. ^ we show, in a log-log scale, w{N,p) as a function of N for many p values, togheter with the fits to the 
data. The uppermost and lowermost lines have slopes —1/3 and —1/2 respectively. Data for p < 0.5 can be perfectly 
fitted by simple power laws (straight lines in Fig. and the resulting z/(p) exponents have been reported in Fig. 
We note that as long as p < 0.25 the exponent turns out to be highly compatible with 3, which is known to be the 
right value for p — 0. Thus we conclude that for p < 1/4 the exponents are those of the p = fixed point. 

For 0.25 < p < 0.5 we find that the exponent takes non-trivial values between 2 and 3. Then one of the following 
two conclusions may hold. Either the transition for p > pq is driven by the p = 1 fixed point and the v exponent is not 
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universal, or more probably any different p value defines a new universality class. This result is very surprising and 
interesting for the possibility that different universality classes are simply the consequence of the random hyper-graph 
topology. 

More complicated is the fitting procedure for p > 0.5. In a recent paper Wilson has shown that in SAT 
problems there are intrinsic statistical fluctuations due to the way one construct the formula. This white noise 
induces fluctuations of order N^^''^ in the SAT probability. If critical fluctuations decay faster than statistical ones 
(i.e. ly < 2), in the limit of large N the latter will dominate and the resulting exponent saturates to = 2. Data for 
p — 0.75 and p = 1 shown in Fig. ^ have a clear upwards bending, which we interpret as a crossover from critical 
(with ly < 2) to statistical — 2) fluctuations. Then we have fitted these two data sets with a sum of two power 
laws, w{N) = AN'^I^ + BN^^I"^ . The goodness of the fits (shown with lines in Fig. ^ confirm the dominance of 
statistical fluctuations for large iV. Moreover we have been able to extract also a very rough estimate of the critical 
exponent v from the subleading term. In Fig. |^ we show with filled squares these values, which turn out to be more 
or less in agreement with a simple extrapolation from p < 0.5 results. 



V. CONCLUSIONS AND PERSPECTIVES 



The exact analysis of a solvable model for the generation of random combinatorial problems has allowed us to show 
that combinatorial phase diagrams can be affected by rare events leading to a mixed SAT/UNSAT phase. The energy 
difference between such SAT and UNSAT instances is non extensive and therefore non detectable by the usual /3 — > oo 
statistical mechanics studies. However, a simple probabilistic argument is sufficient to recover the correct proportion 
of instances. 

Moreover, through the exact location of phase boundaries together with the use of a polynomial global algorithm 
for determining the existence of solutions we have been able to give a precise characterization of the critical exponents 
V depending on the mixing parameter p. The p-dependent behavior conjectured in ref. |^ for the random 2 + p 
SAT case finds here a quantitative confirmation. The mixing parameter dependency also shows that the value of the 
scaling exponents is not completely determined by the nature of the phase transition and that the universality class 
the transtion belongs to is very probably determined by the topology of the random hyper-graph. The model we 
study has also a physical interpretation as a diluted spin glass system. It would be interesting to know whether the 
parameter-dependent behavior of critical exponent plays any role in some physically accessible systems. 

A last remark on the generalization of the present model. With the same analytical techniques presented here, 
one can easily solve a Hamiltonian containing a fraction fk of /c-spin interacting terms for any suitable choice of 
the parameters fk pOf . The case presented in this paper (/2 = 1 — p and /s = p) is the simplest one. There are 
choices which show a phase diagram still more complex with, for example, a continuos phase transition preceded by 
a dynamical one. 
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